Intra-cellular transport of single-headed molecular motors KIF1A 
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Motivated by experiments on single-headed kinesin KIF1A, we develop a model of intra-cellular 
transport by interacting molecular motors. It captures explicitly not only the effects of ATP hy- 
drolysis, but also the ratchet mechanism which drives individual motors. Our model accounts for 
the experimentally observed single molecule properties in the low density limit and also predicts a 
phase diagram that shows the influence of hydrolysis and Langmuir kinetics on the collective spatio- 
temporal organization of the motors. Finally, we provide experimental evidence for the existence of 
domain walls in our in-vitro experiment with fluorescently labeled KIF1A. 
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Intra-cellular transport of a wide variety of cargo in 
eucaryotic cells is made possible by motor proteins, like 
kinesin and dyncin, which move on filamentary tracks 
called microtubules (MT) [HQ. However, often a sin- 
gle MT is used simultaneously by many motors and, in 
such circumstances, the inter-motor interactions cannot 
be ignored. Fundamental understanding of these col- 
lective physical phenomena may also expose the causes 
of motor-related diseases (e.g., Alzheimer's disease) |3( 
thereby helping, possibly, also in their control and cure. 
Some of the most recent theoretical models of interact- 
ing molecular motors HHH0 utilize the similarities be- 
tween molecular motor traffic on MT and vehicular traffic 
on highways Q both of which can be modelled by appro- 
priate extensions of driven diffusive lattice gases [E IH| • 
In those models the motor is represented by a self-driven 
particle and the dynamics of the model is essentially an 
extension of that of the asymmetric simple exclusion pro- 
cesses (ASEP) HE3 

that includes Langmuir-like kinetics 
of adsorption and desorption of the motors. In reality, 
a motor protein is an enzyme whose mechanical move- 
ment is loosely coupled with its biochemical cycle. In 
this letter we conside r sp ecifically the single-headed ki- 
nesin motor, KIFIA [111 Il2l ll-l Il4| ; the movement of 
a single KIFf A motor has been modelled recently with 
a Brownian ratchet mechanism 0, . In contrast to 
the earlier models 0, IE IE °f molecular motor traf- 
fic, which take into account only the mutual interactions 
of the motors, our model explicitly incorporates also the 
Brownian ratchet mechanism of individual KIFIA mo- 
tors, including its biochemical cycle that involves adeno- 
sine triphosphate (ATP) hydrolysis. 

The ASEP-like models successfully explain the occur- 
rence of shocks. But since most of the bio-chemistry is 
captured in these models through a single effective hop- 
ping rate, it is difficult to make direct quantitative com- 
parison with experimental data which depend on such 
chemical processes. In contrast, the model we propose 
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FIG. 1: The biochemical and mechanical states of a single 
KIFIA motor. On the left of the dotted line, KIFIA is bound 
to a fixed position on the MT (state 1), while on the right it 
diffuses along the MT track (state 2). At the transition from 
state 1 to 2, KIFIA detaches from the MT. 



incorporates the essential steps in the biochemical pro- 
cesses of KIFIA as well as their mutual interactions and 
involves parameters that have one-to-one correspondence 
with experimentally controllable quantities. 

The biochemical processes of kinesin-type molecular 
motors can be described by the four states model shown 
in Fig. bare kinesin (K), kinesin bound with 

ATP (KT), kinesin bound with the products of hydrol- 
ysis, i.e., adenosine diphosphate(ADP) and phosphate 
(KDP), and, finally, kinesin bound with ADP(KD) af- 
ter releasing phosphate. Recent experiments |n], 0] re- 
vealed that both K and KT bind to the MT in a stereo- 
typic manner (historically called "strongly bound state" , 
and here we refer to this mechanical state as "state 1"). 
KDP has a very short lifetime and the release of phos- 
phate transiently detaches kinesin from MT ^4|. Then, 
KD re-binds to the MT and executes Brownian motion 
along the track (historically called "weakly bound state" , 
and here referred to as "state 2"). Finally, KD releases 
ADP when it steps forward to the next binding site on the 
MT utilizing a Brownian ratchet mechanism, and thereby 
returns to the state K. 

Model definition. — A single protofilament of MT is 
modelled by a one-dimensional lattice of L sites each of 
which corresponds to one KIFlA-binding site on the MT; 
the lattice spacing is equivalent to 8 nm which is the sep- 
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FIG. 2: A 3-state model for molecular motors moving along 
a MT. denotes an empty site, 1 is K or KT and 2 is KD. 
Transition from 1 to 2, corresponding to hydrolysis, occurs 
within a site whereas movement to the forward or backward 
site occurs only when motor is in state 2. At the minus and 
plus ends the probabilities are different from those in the bulk. 



aration between the successive binding sites on a MT 
0. Each kinesin is represented by a particle with two 
possible internal states labelled by the indices 1 and 2. 
Attachment of a motor to the MT occurs stochastically 
whenever a binding site on the latter is empty. Attach- 
ment and detachment at the two ends of the lattice need 
careful treatment and will be specified below. Thus, each 
of the lattice sites can be in one of three possible allowed 
states (Fig. |3J : empty (denoted by 0) , occupied by a ki- 
nesin in state 1, or occupied by a kinesin in state 2. 

For the dynamical evolution of the system, one of the 
L sites is picked up randomly and updated according to 
the rules given below together with the corresponding 
probabilities (Fig. EJ: 
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2 with Uhdt 

■* 1 with Ugdt 
-> 01 with ujfdt 

20 02 with uj b dt 
02 20 with uj b dt 
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The probabilities of detachment and attachment at the 
two ends of the MT may be different from those at any 
bulk site. We choose a and 5, instead of oj a , as the prob- 
abilities of attachment at the left and right ends, respec- 
tively. Similarly, we take 71 and /3i, instead of u>d, as 
probabilities of detachments at the two ends (Fig. |5J). 
Finally, 72 and /?2, instead of u>b, are the probabilities 
of exit of the motors through the two ends by random 
Brownian movements. 

Let us relate the rate constants ujf,ui s and Lo b with the 
corresponding physical processes in the Brownian ratchet 
mechanism of a single KIF1A motor. Suppose, just like 
models of flashing ratchets [HI ll^ | , the motor "sees" a 
time-dependent effective potential which, over each bio- 
chemical cycle, switches back and forth between (i) a pe- 
riodic but asymmetric sawtooth like form and (ii) a con- 
stant. The rate constant ujh in our model corresponds to 
the rate of the transition of the potential from the form 
(i) to the form (ii) . The transition from (i) to (ii) happens 
soon after ATP hydrolysis, while the transition from (ii) 



to (i) happens when ATP attaches to a bare kinesin [Tl|. 
The rate constant oj b of the motor in state 2 captures 
the Brownian motion of the free particle subjected to 
the flat potential (ii). The rate constants lu s and ujf are 
proportional to the overlaps of the Gaussian probability 
distribution of the free Brownian particle with, respec- 
tively, the original well and the well immediately in front 
of the original well of the sawtooth potential. 

Let us denote the probabilities of finding a KIF1A 
molecule in the states 1 and 2 at the lattice site i at time 
t by the symbols ri and hi, respectively. In mean-field 
approximation the master equations for the dynamics of 
motors in the bulk of the system are given by 



drj 
dt 
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W (l - n - hi) - LU h r, t - LU d r. t + w s /i 4 

- fj - hi), (6) 

-uj s hi + Lu h ri - Wfhi(l - r l+ i - h i+ \) 

~uj b hi(2 - r i+ i - h i+ i - Ti-x - hi-i) 
+oj b (hi-x + /ij + i)(l - n - hi). (7) 



The corresponding equations for the boundaries, which 
depend on the rate constants a, 5, ji and [3i for entry and 
exit (Fig. I2J), are similar and will be presented elsewhere 

From experimental data [ill Il2fl . good estimates for 
the parameters of the suggested model can be obtained. 
Assuming that one timestep corresponds to 1 ms, each 
simulation run had a duration of 1 minute in real time. 
The length of MT is fixed as L = 600. The detach- 
ment rate u>d — 0.0001 ms -1 is found to be indepen- 
dent of the kinesin population. On the other hand, 
u> a = 10 7 C/M-s depends on the concentration C (in 
M) of the kinesin motors. In typical eucaryotic cells 
in-vivo the kinesin concentration can vary between 10 
and 1000 nM. Therefore, the allowed range of uj a is 
0.0001 ms" 1 < uj a < 0.01 ms" 1 . The rate w^ 1 must be 
such that the Brownian diffusion coefficient D in state 
2 is of the order of 40000 nm 2 /s; using the the rela- 
tion u>b ~ _D/(8nm) 2 , we get uib — 0.6 ms -1 . Moreover, 
from the experimental observations that ujf/uj s ~ 3/8 
and w s + ujf ~ 0.2 ms -1 , we get the individual estimates 
u>s ~ 0.145 ms -1 and uif ~ 0.055 ms -1 . The exper- 
imental data on the Michaelis-Menten type kinetics of 
hydrolysis Q suggest that 



0.1 mM 



ATP concentration (in mM) 



ms (8) 



so that the allowed biologically relevant range of lu^ is 
0<u; h < 0.25 ms" 1 . 

Single-molecule properties. — An important test for 
the model is provided by a quantitative comparision of 
the low density properties with empirical results. Sin- 
gle molecule experiments ^lj on KIF1A have established 
that 

(i) v, the mean speed of the kinesins, is about 0.2 nm/ms 
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TABLE I: Predicted transport properties from this model in 
the low-density limit for four different ATP densities, r is 
calculated by averaging the intervals between attachment and 
detachment of each KIF1A. 



if the supply of ATP is sufficient, and that v decreases 
with the lowering of ATP concentration following a 
Michaelis-Menten type relation like (JSJ); 

(ii) D/v ~ 190 nm, irrespective of the ATP concentra- 
tion, where D is the diffusion constant; 

(iii) r, the mean duration of the movement of a kinesin 
on the MT, is more than 5 s, irrespective of the ATP 
concentration. 

The corresponding predictions of our model (see Table UJ) 
for uj a = a = 1.0 x 10~ 6 ms , which allows realization of 
the condition of low density of kinesins, are in excellent 
agreement with the experimental results. 

Collective properties. — Assuming periodic boundary 
conditions, the solutions (rj, hi) = (r, h) of the mean-field 
equations Q in the steady-state are found to be 
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-Q h -Q s - (n s -i)k + Vd 

2K(1 + K) ' 
2K 



(9) 
(10) 



where K = LOd/tOa, = ^h/^f, = tjj s /LJf, and 



d = m s K(\ + k) + (n h + n s + (n s - i)K) . (11) 

The probability of finding an empty binding site on a 
MT is Kr as the stationary solution satisfies the equation 
r + h + Kr = 1 . The steady-state flux of the motors along 
their MT tracks is then given by J = uifh(l — r — h). It is 
interesting to note that in the low ATP concentration 
limit (u>h <C uj s ~ to A of our model, the flux of the 
motors is well approximated by Ji ow = Qesp(^ — p), which 
formally looks like the corresponding expression for the 
totally asymmetric exclusion process, where p is close to 
the Langmuir limit 1/(1 + K) and, 
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as the effective hopping probability [T^. 

Although the system with periodic boundary condi- 
tions is fictitious, the results provide good estimates of 
the density and flux in the corresponding system with 
open boundary conditions, particularly, in the high u a 
regime (Fig. which corresponds to jammed traffic of 
kinesin on MT (see Fig.^J). We also see that, for a given 
uj a , the bulk density of motors in state 2 exceeds that of 
those in state 1 as Uh increases beyond a certain value. 



FIG. 3: The stationary density profiles for cu^ = 0.1 (left) and 
ujh = 0.2 (right) in the case u a = 0.001. The blue and red 
lines correspond to the densities of state 1 and 2, respectively. 
The dashed lines are the mean-field predictions © and (I10P 
for periodic systems with the same parameters. 



Phase diagram. — In contrast to the phase diagrams in 
the a — /3-plane reported by earlier investigators 
we have drawn the phase diagram of our model (Fig.0J) in 
the u a —0Jh plane by carrying out extensive computer sim- 
ulations for realistic parameter values of the model with 
open boundary conditions. The phase diagram shows the 
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FIG. 4: Phase diagram of the model in the w/,-w„ plane, with 
the corresponding values for ATP and KIF1A concentrations 
given in brackets. These quantities are controllable in exper- 
iment. The boundary rates are a — Lo a ,/3i,2 = ^d,7i.2 = 8 = 

0. The position of the immobile shock depends on both ATP 
and KIF1A concentrations. 

strong influence of hydrolysis on the spatial distribution 
of the motors along the MT. For very low ojh no kinesins 
can exist in state 2; the kinesins, all of which are in state 

1, are distributed rather homogeneously over the entire 
system. In this case the only dynamics present is due to 
the Langmuir kinetics. 

Even a small, but finite, rate u>h is sufficient to change 
this scenario. In this case both the density profiles and 
p'j of kinesins in the states 1 and 2 exhibit a shock. As in 
the case of the ASEP-like models with Langmuir kinetics 
[^,|f|, these shocks are localized. In computer simulations 
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FIG. 5: Formation of comet-like accumulation of kinesin at 
the end of MT. Fluorescently labeled KIF1A (red) was intro- 
duced to MT (green) at 10 pM (top), 100 pM (middle) and 
1000 pM (bottom) concentrations along with 2 mM ATP. The 
length of the white bar is 2p.m. 

we have observed that the shocks in density profiles of 
kincsins in the states 1 and 2 always appear at the same 
position. Note that if the individual density profiles pj 
and pj exhibited shocks at two different locations, two 
shocks would appear in the total density profile pj = 
p]+Pj violating the usual arguments 01 that ASEP-type 
models exhibit exactly one shock. Moreover, we have 
found that the position of the immobile shock depends 
on the concentration of the motors as well as that of ATP; 
the shock moves towards the minus end of the MT with 
the increase of the concentration of kinesin or ATP or 
both (Fig. HI). 

Finally, we present direct experimental evidence that 



support of the formation of the shock. The "comet-like 
structure" , shown in the middle of Fig. [5] is the collec- 
tive pattern formed by the red fluorescent labelled ki- 
nesins where a domain wall separates the low-density re- 
gion from the high-density region. The position of the 
domain wall depends on both ATP and KIF1A concen- 
trations. Moreover, as we increase the concentration of 
KIF1A, the transition from the regime of free flow of ki- 
nesins to the formation of the shock is observed(top and 
middle in Fig. [5]). Furthermore, we observe jammed traf- 
fic of kinesins at sufficiently high concentration (bottom 
in Fig. 0). The position of the shock in our simulation 
agrees well with the location of the domain wall in the 
comet-like structure observed in experiments^]. 

In this letter we have developed a stochastic model for 
the collective intra-cellular transport by KIF1A motors, 
by taking into account the biochemical cycle of individual 
motors involving ATP hydrolysis and their mutual steric 
interactions. We have been able to identify the biolog- 
ically relevant ranges of values of all the model param- 
eters from the empirical data. In contrast to some ear- 
lier oversimplified models, the predictions of our model 
are in good quantitative agreement with the correspond- 
ing experimental data. Moreover, we have mapped the 
phase diagram of the model in a plane spanned by the 
concentrations of ATP and KIF1A, both of which are 
experimentally controllable quantities. Finally, we have 
reported the experimental observation of a comet-like col- 
lective pattern formed by the kinesin motors KIF1A and 
identified the domain wall in the pattern with the shock 
predicted by our model. 
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